Multiphoton electron detachment by a superposition of static and ac fields 
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A theory of electron detachment from atoms or negative ions by a superposition of a static 
and a laser (or, more generally, ac) fields, parallel to each other, is developed in the case when the 
photon energy and the field amplitude are much less than the detachment energy and intra-ion field, 
respectively. Simple analytical results together with their validity conditions are obtained. Several 
qualitatively different regimes of detachment have been identified. Applications of the present theory 
to electron emission from metal and semiconductor surfaces and from photosynthetic bio-complexes 
are discussed. 
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I. INTRODUCTION 

The photoefFect by a strong nonresonant field in the 
presence of a static electric field is a problem of a sig- 
nificant scientific interest. This problem can arise in dif- 
ferent situations. A static field can accelerate ionization 
of atoms and electron detachment from negative ions in 
the presence of a nonresonant laser field It is also 
of interest to consider effects of an ac field on static-field 
induced electron emission from metal and semiconductor 
surfaces and from photosynthetic bio-complexes. Such 
experiments can be done, e.g., by a modification of the 
techniques of scanning tunneling microscopy 0-111 • 

Such phenomena are not completely understood yet, in 
spite of the fact they have been investigated in a number 
of studies. Ionization of atoms by a strong nonresonant 
laser field in the presence of a static field was first consid- 
ered in Ref. [T] , in which the ionization rate was obtained 
with exponential accuracy (i.e., without preexponential 
factor), in the limits of tunneling and multiphoton ioniza- 
tion. The general exponential-accuracy solution, which 
includes the above limits as special cases, was obtained 
in Ref. @. Attempts to obtain a complete solution for 
electron detachment from negative ions were made in 
Refs. 0-11 ■ However, the resuhs obtained in Refs. 0-11] 
involve infinite series and are quite intricate. Moreover, 
Ref. is concerned with the relatively simple, but prac- 
tically not very interesting, case (a circularly polarized 
light perpendicular to the static field). The analysis in 
Ref. jQj was limited to several simple cases where only one 
term of an infinite series was taken into account. There 
are also direct numerical solutions of the problem p^ . 
However these laborious calculations cannot provide the 
dependence of the ionization rate on the parameters for 
the entire parameter space. Thus, at present there does 
not exist a comprehensive analysis of multiphoton de- 
tachment in the presence of a laser (or, more generally, 
ac) field for the whole range of the parameters of the 
problem 

There has also been recent research on one- ITll - tTsj or 
few-photon (l3 - tl^ ionization in the presence of a static 



field. In particular, a ripplelike structure in the pho- 
todetachment cross section was observed [ll| , which was 
explained by interference of electron waves reflected from 
the barrier due to the static field. 

In this paper, we consider photodetachment from 
atoms or negative ions produced by a static and a non- 
resonant ac fields. We focus on the case of a linearly po- 
larized ac field parallel to the static field, but the present 
treatment can be extended to the general case of elliptic 
polarization and arbitrary field orientations. We provide 
analytical results for the differential and total rates of 
detachment. As an application of the present theory, we 
discuss electron emission from metals and semiconduc- 
tors and from photosynthetic bio-complexes in the pres- 
ence of a static and an ac fields. 

The present theory is an extension of the Keldysh the- 
ory of ionization and electron detachment by a strong 
laser field 17|. The Keldysh theory provides the ion- 
ization and detachment rates and the energy spectrum 
of the emitted electrons in simple analytical forms. The 
nonlinear effects of strong field are taken into account 
with the help of Volkov functions, which describe exactly 
the electron-field interaction in the final electron state. 
The Volkov function for an electron in a harmonic field 
has, at least, two important advantages. First, it has 
a simple analytical form that is convenient for calcula- 
tions. Second, the probability distribution of the ioniza- 
tion products in the Volkov-function basis directly pro- 
vides the energy spectrum of the products. This is due 
to the fact that in the limit of vanishing laser field the 
Volkov function tends to an energy eigenstate of a free 
particle, i.e., a plane wave. 

The standard Volkov functions, obtained for the case 
of a harmonic field, cannot be used in the present case. 
Hence an important step in the solution of the present 
problem is to find suitable functions providing the basis 
for the expansion of the final state of the electron. An 
obvious candidate for such functions are the Volkov func- 
tions, describing an electron in a static and an ac fields 
[7|. These Volkov functions have a rather simple form. 
However, they have the disadvantage that in the limit of 
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the vanishing ac field, they do not reduce to eigenstates 
of the Hamiltonian for a particle in a uniform static field. 

There is also an alternative basis, consisting of func- 
tions which have the important property that they be- 
come eigenstates of the above Hamiltonian in the limit 
of vanishing ac field Q . These functions are expressed in 
terms of the Airy function and, as a result, unfortunately, 
are less convenient for calculations than Volkov functions. 
Therefore, we use a combined approach, in which initially 
we calculate the final state of the electron in the basis of 
the Volkov functions, and then we transform this state 
into the second basis. This approach yields simple an- 
alytical expressions for the rate of electron detachment 
and the energy spectrum of the emitted electrons. These 
results allow us to obtain a comprehensive picture for 
the behavior of the detachment rate for a broad range of 
the parameters of the problem. Our theory agrees with 
Keldysh theory [l3| in the limit of the vanishing static 
field. 

The interaction of matter with electromagnetic fields 
can be described using different gauges. Correspond- 
ingly, in analytical treatments of multiphoton and above- 
threshold ionization [l^ various gauges have been used. 
Most commonly, the so called velocity and length gauges 
have been used. Although the exact solutions are gauge- 
independent, in approximate analytical theories, different 
gauges yield different results. There are sound theoretical 
considerations [IQ] as well as empirical evidence [13, Hlj 
indicating that the length gauge has an advantage over 
other gauges. In the present paper, following a number 
of previous studies, such as, e.g., Ref. we use the 
length gauge. 

In Sec. |ll]the general theory of bound-free transitions 
in a time-dependent field is reviewed. In Sec. IIIII we 
introduce two sets of final-state wave functions: the gen- 
eralized Volkov functions and functions that are a direct 
extension of the stationary states for a particle in a static 
uniform field. In Sec. IIVI we consider the limit of a very 
slowly varying field. Section |V] presents a full solution 
of the problem. In contrast to previous studies, we use 
the stationary-phase method, which provides a relatively 
simple analytical solution. In Sec. |VT]we describe sev- 
eral important regimes of detachment, which follow from 
the above solution. Our solutions are much reacher than 
solutions in the presence of only one of the two fields, 
and a number of novel features are revealed. In Sec. IVIII 
we discuss application of the present theory to electron 
emission from solid surfaces. Section rVIIII provides con- 
cluding remarks. In Appendix [Xj we obtain the expec- 
tation value of the energy in various states describing an 
electron affected simultaneously or separately by static 
and time-dependent fields. In Appendix |B1 we calculate 
a useful integral involving the Airy function. 



II. QUANTUM TRANSITIONS DUE TO A 
NONSTATIONARY FIELD 

In this section we summarize our general approach, 
which is similar to that developed by Keldysh, Faisal, 
and Reiss fl7l. [2^ [23j for a description of above-threshold 
ionization [18]. The Hamiltonian of a quantum system 
(e.g., an atom or ion) in a field is 

H{t)^ Ha + Vit), (2.1) 

where Ha is the atom Hamiltonian and V{t) describes 
the atom-field interaction, which induces transitions to 
the continuum from the initial bound atomic state |(/>o), 
an eigenstate of Ha with the energy Eq. We assume that 
only one electron interacts effectively with the field. The 
field is assumed to be sufficiently weak, so that it does not 
modify appreciably the initial state and its energy. More- 
over, the field alternates periodically with the frequency 
Lo and is nonresonant. In view of the above assumptions, 
we can neglect transitions to other discrete states and 
cast the wave function in the form, 

\^f{t)) = m\M+J2^>^(^^\Mt)), (2.2) 

A 

where the functions 1-0^(0) ioim an orthonormal basis of 
the continuum and I3{t) and /3x{t) are the amplitudes of 
the initial and final states, respectively. 

For t ^ 27r/a;, the initial-level population is given by 

\m\' « e-'^*. (2.3) 

where W is the rate of bound-free transitions. The rate 
W is obtained most easily by considering short times t ^ 
VF-i, when Eq. (gS]) becomes |/3(i)p « 1 - PFi. Then 
the normalization condition for the wave function (j2.2p 
implies that 

^ = ,1™ 7El/^A(i)p. (2.4) 

A 

To obtain I3\{t) in Eq. (|2.4p . we use the time-dependent 
perturbation theory. Namely, we note that for short 
times, in the first approximation, 

m~exp(^-^Eoty (2.5) 

Assume that the functions \tp\{t)) are the solutions of the 
Schrodinger equation 

ih^^ = HFit)\,Px). (2.6) 

Here Hp{t) is the Hamiltonian for the free (unbound) 
electron in the field. Then one can show with the help of 
(1231) that [H m 

m) = - J / dt'{Mt')\V{t')\cf>o) exp (~^Eot'^ ■ 

(2.7) 
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In this paper, we use the common approximation [13, [H, 
[23} in which the interaction of the detached electron with 
the field is taken into account, but the interaction with 
the atom is neglected. Then 

HF{t)^^ + V{t), (2.8) 

where m is the electron mass and P is the momentum 
operator. 

Note that the basis {\i^x{t))} satisfying condition p.6|) 
is not unique, being defined with the accuracy to an ar- 
bitrary time-independent unitary transformation. For a 
system with a time-independent Hamiltonian, there is a 
preferable basis consisting of the eigenstates of the Hamil- 
tonian. However, in the time-dependent case, such a 
preferable set of solutions of the Schrodinger equation 
is absent (the eigenstates of the Hamiltonian are not so- 
lutions of the Schrodinger equation) , and therefore more 
than one basis {\'>P\{t))} may appear convenient, depend- 
ing on the situation. In particular, in the present formal- 
ism, two types of states of the electron in a field prove to 
be useful. These states are described in the next section. 

III. FINAL-STATE WAVE FUNCTIONS 

Here the field-atom coupling is described in the dipole 
approximation and in the length gauge; correspondingly, 
the interaction term in Eqs. ()2.1|) and ()2.8|) is 

Vir,t) = -f{t)-f. (3.1) 

Here f is the radius- vector of the electron initially bound 
to the atom, and ^{t) is the force exerted on the electron 
due to a uniform electric field J^(t)/eo, where Cq is the 
(negative) electron charge, eg = — |eo|. 

In this paper, the electric field J-{t)/eo is assumed to 
be a superposition of a dc term £ /cq and a harmonic field 
with amplitude F/cq, 

f{t) = £ + Fsmujt. (3.2) 
Below we use the units in which 

h = m = l. (3.3) 

A. Volkov functions 

1 . General formulas 

Equation (|l!6)) with the account of Eqs. (IZl)) and ((3?T| . 
where J-{t) is an arbitrary function, has solutions in the 
form of plane waves, 



where 

P{t)^p-^A{t), (3.5) 
Here c is the speed of light and A{t) is defined by 

A^~cf{t)/eo. (3.6) 

It is not difficult to check that the function p.4|) indeed 
satisfies Eq. (|2.6p . The functions of the form given in 
Eq. ((Ol) are called the Volkov function [13, Hi]. 

The function ijjp{r, t) is an eigenstate of the momentum 
operator with the eigenvalue Pit) in Eq. p.Sp . Note that 
in the length gauge assumed here, the vector potential 
vanishes, and as a result the momentum coincides with 
the kinetic momentum p.Sp . The vector p in Eq. (|3.5p is 
a constant of motion. 

The Volkov functions satisfy the orthonormality con- 
dition 

{^p{t)\i^p(t))^{2nf6{p~^) (3.7) 
and the completeness condition 

/ ■^^'W^t)Mr.t) = <5(f- r-). (3.8) 

Therefore the Volkov functions for a given F{t) form a 
basis of the Hilbert state of a free electron. 



2. Volkov functions for a sum of static and ac fields 

The general Volkov functions p.4p are applicable, in 
particular, for the field given in Eq. p.2p 25]. Now in 
Eq. dSH), P{t) is described by Eq. with 

en - F 

— A(t)^-£t+—cosujt. (3.9) 

c UJ 

In the present paper, we focus on the important case, 
in which the ac electric field is linearly polarized in the 
direction of the constant electric field. On choosing the 
latter direction as the z axis, the interaction term p.ip 
involved in Eqs. (j2.ip and (|2.8p becomes 

V{z,t) = -T{t)z, (3.10) 

where the force J-{t), as implied by Eq. (|3.2p . is given by 

F{t) = E + Fsinujt {£,F>0). (3.11) 

Here, without loss of generality, we assume that £ and F 
are positive. Now 

A{t)^ (0,0, A{t)), (3.12) 

where 

en F 

—A{t) = ~£t+—cosu}t, (3.13) 

c CJ 
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and, correspondingly, Eq. (13.51) becomes 



with 



F 

Pz{t) = Pz + St cosujt 

UJ 



(3.14) 



(3.15) 



Thus, the transverse momentum q = (j)x,Py) is con- 
served, whereas the z-component of the momentum 
varies with time, the quantity pz being a constant of mo- 
tion. 

Then Eq. (|3.4p becomes 

' p^ 



^p{f, t) = exp <^ i 



^ ,2 



£p 

2 

EF 



£2 



q-p + Pz{t)z- |^£^ + C/p)t 

FPz . . ^ 

— — sm ojt sm 2ujt 



-t sinujt 



£F 



■ cos Lot 



(3.16) 



where p — {x, y) and Up is the ponderomotive energy, 
i.e., the average kinetic energy of the electron oscillating 
in the harmonic field f29l. 



4cj2 



F^ 
4ma;2 



(3.17) 



the latter expression being written in the usual units. In 
the derivation of Eq. (|3.16p , we omitted a constant phase. 

The Volkov functions p.l6p have a relatively simple 
form and hence can be convenient for analytical calcu- 
lations. However, in the next subsection we show that 
the Volkov functions do not reduce to eigenstates of the 
energy operator in the limit of the vanishing ac field, and 
hence the Volkov functions are not suitable for obtaining 
the energy spectrum of the detached electrons. There- 
fore, in Sec. IIII CI we introduce another set of solutions 
of Eq. (HH) with the account of Eqs. (g^, (|XTU)) . and 
p.ll[) . as a generalization of the stationary states for a 
particle in a uniform field. 



B. Static field: Stationary states and Volkov 
functions 

In this subsection, we consider two types of the electron 
states for the case of a static, uniform field. 



1. Stationary states 

The stationary states of a particle in a static, uniform 
field have the form [2^ 

^Ed^ = ^EA^)e''-', (3.18) 
where E is the energy of the state and 



2i/3 

i^EAz) = ^T7^Ai(-M) 



(3.19) 



is the energy eigenstate of one-dimensional motion of a 
particle in a uniform field. Here 



u= {2£y^^{z + Ez/£), Ez=E-q^/2 (3.20) 



and 



m) = ]^ dvcos(^^+^v^ (3.21) 

is the Airy function [l^, . Note the asymptotic expres- 
sions for the Airy function [27j . 



Ai(0 



1 



expl--C3/2), ^»1 (3.22a) 



Ai(C) = (^1^1'/' + ' » 1' (3-22b) 

whereas Ai(0) w 0.36. Thus, the Airy function decays 
exponentially for positive values of the argument and os- 
cillates for negative values of the argument. 



2. Volkov functions for a static field 

Along with the standard basis of the stationary func- 
tions (I3.18P described in textbooks, we consider an alter- 
native basis of the static-field Volkov functions, following 
from Eq. for = 0, 



V'p(r, t) = exp {i[q- p + {pz + £t)z 



'_t-^t^-—t 
6 



(3.23) 



The two bases are related by a unitary transformation of 
the form 



-0j?(r, t) 



dE hq{pz,E)i)Eq{r)e 



iEt 



(3.24) 



The function hq {pz , E) here can be obtained as follows. 
Equation (|3.24p has the form of a Fourier transform and 
hence can be inverted: 



hq{pz,E)lljEq{r) 



1 

2^ 



dt ipp(r, t)e 



iEt 



(3.25) 



This integral can be transformed to the form proportional 
to Eq. (I3.18P with the coefficient 



hq{pz,E) ^ -^^^P 



£ \2 6 



(3.26) 



It is easy to see that the function in Eq. (I3.26|) satisfies 
the relations 



dp. 



2tt 



'h;ipz,E')hq(pz,E) = S{E-E'), (3.27) 



dEhqipz,E)h;{p'z,E)^2nS{pz -p'z), (3.28) 
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which ensure that the transformation in Eq. (|3.24p is uni- 
tary. Using the unitarity of hq(pz,E), we can invert the 
transformation in Eq. p.24|) . as foUows. Multiplying the 
both sides of Eq. (IXM)) by h*{p^,E')/{2n), integrating 
over pz, and using Eq. (I3.27p . we obtain 



-iEt 



dpz 
2tt 



The Volkov function in Eq. p.23p has a relatively sim- 
ple form, but it is not an eigenfunction of the energy oper- 
ator. In fact, the energy distribution in a Volkov state is 
infinitely broad. Indeed, as follows from Eq. (|3.24l) . this 
distribution is proportional to \hq{pz, E)\'^ and hence is 
uniform [cf. Eq. p.26p ]. Note, however, that, as shown 
in Appendix lA 11 the expectation value of the energy in 
the state is finite and equals p'^/2 [see Eq. (PT6)) ]. 



yields 
i^Eq{r, t) 



21/3 



exp < I 



q- p 



■ sin 2ijjt — Et 



F 

z I — cos (jjt 
M{-u{t)), (3.32) 



h*q{p,,E)^:p{r,t). (3.29) where 



uit) = (2£)i/3 + 



F 



■ s'm ujt 



(3.33) 



As shown in Appendix I A 21 [Eg. (jA24p ]. the mean energy 
in the state \ipEq{t)) averaged over the oscillations of the 
electron is equal to E. 

One can invert the unitary transformation p.30p with 
the help of Eq. (1X28)) to obtain 



\M-t))= dEhqip,,E-Up)\i^Eq{t)). (3.34) 



Alternative basis of exact solutions for an 
electron in static and ac fields 



As mentioned above, the Volkov functions are exact 
solutions of Eq. (|2.6p . and they constitute a basis for the 
Hilbert space of a free electron. This basis is not unique — 
actually, any transformation of the Volkov functions with 
the help of a time-independent unitary operator produces 
a basis consisting of solutions of Eq. (|2.6p . Here we in- 
troduce one such basis, consisting of functions which, in 
contrast to the Volkov functions, reduce to stationary 
states in the limit of the vanishing ac field. 

In the case of a time-dependent field p. Ill) , we consider 
a direct extension of the unitary transformation p.29p . 



li^Eqit)) 



dpz 
2tt 



h:{p,,E-Up)\Mt)), (3-30) 



where |V'p(0) the Volkov function for a particle in the 
field (PTT|) [see Eq. ([XTC]) ]. Equation P^UI defines a 
new basis consisting of solutions of Eq. (|2.6I) with the 
account of Eqs. (P?51) . (IXTU|) . and (I5TT|) . in addition to 
the basis of the Volkov functions in Eq. (|3.16p . Taking 
into account Eqs. (|3.7p and (I3.27p . one can show that the 
fimctions p.30p satisfy the normalization condition 

{i'EqimE'qAt)) = {2TTfS{q-q')6{E-E'). (3.31) 



In Eq. (|3.30p , the energy is shifted by the value of the 
ponderomotive energy, Up. This corresponds to the fact 
that the averaged motion of a particle affected simulta- 
neously by a time-independent and a rapidly oscillating 
fields occurs as if in addition to the time-independent 
field there is another time-independent field with the po- 
tential energy equal to the ponderomotive energy ,29] . 

To obtain the state I'^I^Eqit)) in the coordinate repre- 
sentation, in Eq. p. 301) we insert Eqs. p. 261) and replace 
[■(/'p'(t)) by the function in Eq. (I3.16P : then the integration 



IV. IONIZATION AND DETACHMENT IN THE 
ADIABATIC LIMIT 

A. Units of measurement 

Henceforth, we use dimensionless units [s^ defined by 
relations p.3p and by the equality 



Eo = -1/2. 



(4.1) 



When the initial state is the ground state of a hydrogen 
atom, the present units coincide with the familiar atomic 
units where h = m = Cq = 1. However, generally the 
present units differ from the atomic units. 

The transition in any expression from the present di- 
mensionless units to the conventional units is performed 
by the substitution 



(4.2) 



for each quantity X entering the expression, where Xa 
is the unit of X defined by Eqs. p^]) and (l4ll . For 
example, the units of frequency, momentum and force 
are 



UJa = 2\Eo\/h, 
Pa = v^2m|£;o| 



Fa = LOaPa = V^i^lEoD^'/yn 



Ionization in a static field 



(4.3a) 
(4.3b) 
(4.3c) 



Before developing a general theory, we consider the 
simple special case of a field with a very low frequency 
(the adiabatic limit). We start with ionization or detach- 
ment of an atom (or ion) by a static field T/bq. In this 
case the ionization rate is given by 



C(F)e 



-2/(3|.F|) 



(4.4) 
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where C{T) is the preexponential factor which depends 
on the specific details of the problem. 

For instance, for field-induced ionization of a hydrogen 
atom in the ground state [3l| 



C{T)=A/\T\, 



(4.5) 



whereas for electron detachment from an s-state in a po- 
tential well due to a short-range potential, such as, e.g., 
the potential in a negative ion, one has [30| 



C{J^)^ttA^\J^\. 



(4.6) 



Here A is the real coefficient in the asymptotic expression 
for the electron wave function. 



(4.7) 



where rg is the radius of the potential (tq ^ 1 ). I n par- 
ticular, for the case of a weakly bound s-state [32| 



A=(27r)-i/2 (ro«l). 



Static and ac fields: Adiabatic limit 



(4.8) 



If the field slowly varies with time, Ec^. ()4.4|) is as- 
sumed to be still valid (the adiabatic limit). In the case of 
the periodic field J^{t) given in Eq. (|3.1ip . the transition 
dynamics averaged over the field period is exponential 
with the rate obtained by averaging Wst over the period, 

W:^—f dwC(£-|-Fsinw)e-2/(3l^^+^''^i"''l). (4.9) 

This holds if the field frequency is much greater than the 
average transition rate. 



(4.10) 



Let us calculate the integral in Eq. ()4.9p . The expo- 
nential in the integrand of Eq. (|4.9p has the absolute 
maximum aX v — 7r/2, where the field value is maximal, 
which occurs at times {tt/uj){21 + 1/2) {1 = 0,1,...). This 
maximum is the only local maximum for £ > F , whereas 
for £ < F, there is also a maximum at v = —ixjl, which 
is lower than the maximum at u = 7r/2. In the case 
considered here when £ affects significantly bound-free 
transitions, the contribution of the second maximum is 
relatively small and hence will be neglected. Using the 
saddle-point method, we obtain 

W = \\[^{£ + F)C{£ + F)e-2/[3(£+^^)]. (4.11) 
2 V T^F 

In particular, for photodetachment we obtain, with the 
account of Eq. 



W ■ 




(£ + F)2e-2/[3(£+^^)l. 



(4.12) 



For validity of Eq. (|4.11l) , the width of the peak in the 
integrand of Eq. (j4.9p , given by the inverse square root of 
the magnitude of the second derivative of the exponent 
in (14. 9p at w = 7r/2, should be much less than one. This 
holds under the following conditions: 



£,F 



(4.13) 
(4.14) 



Inequalities (j4.13l) coincide with the validity condition 
of (|4.4I) . which states that the external field should be 
much weaker than the intraatomic field. In this case the 
field-induced barrier is so broad that the transition rate 
is exponentially small. 

An analysis shows that for F ^ £'^ , the rate in 
Eq. (|4.9p is practically independent of the ac field and 
is approximately given by Eq. (j4.4l) with T = £, whereas 
for 



F>£'' 



(4.15) 



the rate W significantly depends on the ac field. Thus, 
the condition (|4.14p (obtained previously in [i|), in 
essence, coincides with the region in which the ac field 
significantly affects the transition rate. Similarly, when 
£ is smaller than F , the static field significantly affects 
the transition rate when [cf. Eq. (I4.12p ] 



£>F'' 



(4.16) 



Below we will recover Eq. (I4.12p as a special case of 
our theory for sufficiently small w. 



V. DETACHMENT BY THE SUM OF DC AND 
AC FIELDS 

Here we apply the general theory of Sec. [H] to de- 
tachment of negative ions by the field p. lip . Using the 
Volkov functions, the wave function (12. 2p in the coordi- 
nate representation takes the form 



^'(r,t) =/3(t)0o(^) + 



dp 
(2^ 



Pp-{t)^Pp{r,t). (5.1) 



The amplitude of the transition from the initial state to 
a Volkov function follows from Eq. (12.71) with the account 
of Eqs. ([SJO]) and (gl]), 

Pj^it) dt' J dr(t>o{r)F{t')zi!t{r,ty'/^. (5.2) 



On changing the order of integration and using Eqs. p.4p , 
dsn, and (IXTCj) . Eq. dOl) becomes 

j3j;{t)=i J drMr)ze~"^pJ^ dt' F{t')e'f^''\ (5.3) 
where 

f{t) -\ + \f/^ W - + (5-4) 
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A. Evaluation of the time integral in Eq. (|5.3 

1. Saddle points 

We proceed to evaluate the time integral in Eq. ([57 
by the stationary-phase (or saddle-point) method. The 
saddle points are the values of t which satisfy the equa- 
tion 



fit) 



PKt) 



F{t)z- 



2 



0. 



(5.5) 



Here we used Eq, 
Eqs. (1331) and dill), 



and the relation, following from 



(5.6) 



Equation (|5.6p has the form of Newton's second law; this 
fact is a consequence of the Ehrenfest theorem. 

Equation (15. 5p is rather complicated. Fortunately, as 
shown below, the coordinate z and the momentum q are 
small in Eq. (j5.5l) . and hence in the zero-order approx- 
imation we can neglect the last two terms in the right- 
hand side of Eq. (|5.5p . yielding P^{t) = — 1 or, in view of 
Eq. 



F 

Pz{t) = pz + £t cos = is, 



(5.7) 



where s = ±1. Equation (|5.7p implies that the saddle 
points are complex. 

To solve Eq. (15. 7p . we write t as the sum of real and 
imaginary parts, t = tr + iti, and separate the complex 
Eq. (|5.7p into two real equations 



F 

Pz + £tr COS LOtr COSh LOti 

(jj 



£ti H siwuitr smhujti — s. 



0, 



(5.8a) 
(5.8b) 



We assume that the quantity Pz + £tr in Eq. (|5.8ap is 
small (the validity conditions for this assumptions are 
obtained below); hence it can be neglected in the zero- 
order approximation. As a result, Eq. (|5.8al) yields the 
solutions tr = Tfc, where 



rfc = -lfc+i) (fc = 0,l,...) (5.9) 



are the moments of the field extrema. Inserting tr ~ 
into Eq. (j5.8bp yields ti — r], where r/ is the solution of 
the equation 



1] = - - (-1) — smhwT]. 
b Lot 



(5.10) 



Thus, we obtain that in the zero-order approximation the 
saddle points t^ are given by 



t 



(0) 



Tfc -I- ir]. 



(5.11) 



As follows from Eq. (15. 5p . in the zero-order approxi- 
mation 



/(0)(40)) ^ Pz{tf )Pz{tfy) = ^sHtT) = ^sb,. (5.12) 

Here in the second equality we used Eqs. (15. 6p and (|5.7p . 
whereas, in view of Eqs. (13. lip . (|5.1ip . and ([S? 



where 



Ci = Fcoshw?7 = y/WTuSHs~~£^^ 



(5.13) 



(5.14) 



The second equality in Eq. (I5.14p is implied by Eq. (jS.lOp . 

The quantity Pz + £tr, neglected in Eq. (|5.8ap . equals 
to a first approximation 



Pk = Pz +£Tk. 



(5.15) 

Equations (|5.7p and (|5.9p imply that pk = Pz{Tk), i.e., 
the quantity pk is the value of the z-component of the 
momentum at the moment Tk for an electron detached in 
the vicinity of the A:th saddle point. The quantities z, q, 
and Pk are the small parameters of the expansions below. 

To obtain corrections to Eq. ([J" 
of Eq. (|5.5p in the form 



(liOT]), 



we write the roots 



tk = t)"' + t 



(0) , ,(1) 



fe 



(5.16) 



(n) 

where tj^ is the nth-order contribution into tk- Then 
we replace t by Eq. (j5.16p in Eq. ()5.5p with the account 
of Eqs. p. lip and p.l5p . expand the trigonometric func- 
tions into the Taylor series in a neighborhood of tf^ , and 
equate each sum of the terms of the same order to zero. 
In particular. 



This yields 4"^- 



tl 



(2) 



-isz 

-r2 



-Pk/h, 
Si ,\ Sipl.sq' 



where 



Si = {-!)'' Fojsmhuj7] = oj^{s - £t]), 



(5.17a) 
(5.17b) 

(5.18) 



the latter equality following from Eq. (|5.10p . 

Furthermore, we calculate f{tk)- To this end, we note 
that the explicit form for /(t) in Eq. (15.41) is given by the 
expression in the square brackets in Eq. p. 161) . where 
q • p should be replaced by t/2. (In this expression a con- 
stant term is neglected, which provides an unimportant 
constant phase.) Inserting in this expression Eq. (|5.16p 
instead of t and performing the Taylor expansion of the 
resulting expression about t^k'^ , we obtain the expansion 
of/(tfe), 



f{tk) = f^''\tk) + f'^'Htk) + 



(5.19) 



where /'■"■'(tfc) is the contribution of the 7ith order of 
magnitude. In particular. 



f^'^itk) 



+ fc/p + ^ )rfe + zCo, (5.20a) 
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f(^){tk)=isbk=- + iTk+iCk)i^ + iTk+iv)^, (5.20c) 



Pk, (5.20b) 



pI 



■ \1 



where 



(5.21) 



Ck = V- s/h- 



(5.22) 



The above expansions (|5.16l) and (|5.19p are vahd, if the 
quantities z, q, and pk are sufficiently small. As implied 
by the imaginary part of Eq. (j5.20cp [cf. the last integral 
in (|5.3p ]. the coordinate and momentum of the electron 
at the moment have the bounds 

\z\<{sbk)-'/', k|<c-^/^ \q\<V-'^''- (5.23) 

Thus, the condition that q, pk, and z have bounded distri- 
butions (which is a weaker requirement than the small- 
ness of z, q, and pk) is equivalent to the simultaneous 
inequalities 



(a) sbk > 0, (b) 7] > s/bk, 



(5.24) 



and 77 > 0. Condition (a) in Eq. ()5.24p ensures that the 
integrand of the time integral in (j5.3p has a maximum at 
tk [cf. Eq. (|5.12p ] , implying that all roots of Eq. (15.51) are 
saddle points, whereas the condition i] > 0, which is a 
consequence of Eq. (|5.24p . requires that all saddle points 
are located in the upper half space of the complex t plane 
[cf. Eq. dsn])]. 

The conditions in Eq. ()5.24p imply the value of s, as 
follows. Inequality (b) in Eq. ()5.24p can be recast as 77 > 
l/(s6fc), since = 1, or sbkij > 1, in view of condition 
(a) in Eq. ((S^M)) . Hence, we can write 



sbkTj — 1 = sSrj — 1 + {—l)'^srjF coshcjTy 



(— l)'^s77Fcosh(u;7/ ( 1 



tanhwr; 
Lorj 



>0, (5.25) 



where we used Eqs. (j5.13p and (|5.14p in the first equality 
and Eq. (I5.10p in the second equality. Since the last 
factor in the parentheses in Eq. (|5.25l) is positive and 
> 0, the inequality in Eq. (|5.25p yields 



(5.26) 



For the moments Tk with even (odd) fc, the ac field is 
directed along (opposite to) the static field. Correspond- 
ingly, the cases with even and odd k play different roles in 
the process of detachment. For even fc, Eq. (I5.13P implies 
that bk > 0, and condition (a) in Eq. (I5.24p always holds. 
However, for odd fc, when F < £ and uj is sufficiently 
low, condition (a) in Eq. (|5.24p is violated, which means 
that there are no saddle points in such cases. [Indeed, 



then Eq. (I5.13P implies that &fc > 0, whereas s = — 1, in 
view of Eq. (|5.26p . yielding a violation of condition (a) 
in Eq. (I5.24p .] This fact can be easily understood, if one 
takes into account that for F < £ the moments Tk with 
odd k correspond to the field minima^ and hence such 
moments cannot contribute appreciably to detachment, 
at least, in the adiabatic (very low oj) limit (cf. Sec. lIV Cl) . 
Though the case of odd k with F > £ involves local max- 
ima of the field magnitude equal to F ~ £, these maxima 
are less than the maxima equal to F + £, obtained for 
even k. Hence, the case of odd k produces a negligi- 
ble contribution to the detachment rate, unless £ is very 
small [33| . 

Therefore, below we restrict the consideration to the 
saddle points corresponding to the field global maximum, 
i.e., to the case oi k — 21 {I = 0,1, . . .). Then s = 1 [see 
Eq. (j5.26p ]. and the time integral in Eq. (|5.3p is, in view 
of Eq. (|5:T2|) . 

/ dt'J-(t')e^^(*')«^^(4?')e*^(*-) 

X / di'e-''"^*'-*-)'/^ = y2^^e^/(*-\ (5.27) 

where N = [ujt/2TT — 1/4], [. . .] denotes the integer part, 
and 



bo^b2i=£ + Cu Cl = ^F^+uj^{l-£7j)^. (5.28) 



Note that now in Eq. ()5.22p Ck = C21 = cq, where, in view 
of (lOSD . 

co = r?-l/&o. (5.29) 



2. Asymptotic expansions for some parameters 

The present theory involves a number of auxiliary pa- 
rameters, which enter the analytical results. It is useful 
to have simple formulas, providing approximations or es- 
timates of the parameters of the present theory. Here we 
provide asymptotic expansions for some parameters. 

We begin with 77, which is a solution of the transcen- 
dental equation (I5.10p . In the case of even k, Eq. (|5.10l) . 
in view of Eq. (j5.26p . becomes 



£rj = 1 sinhcjr;. 

UJ 



(5.30) 



This equation has a unique solution, which decreases with 
Lj. It is possible to obtain asymptotic expansions for 77 in 
different limits. 

In particular, for small u, one can expand sinhaj77 in 
Eq. ()5.30p in powers of lot], which allows us to obtain the 
expansion. 



1 



Fuj' 



F{£_-^FyJ 

£ + F~ Q{£ + FY^ 120(5 + F)7 



-fO(w^). (5.31) 
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This series is rapidly converging, and hence can be trun- 
cated at any term to a good approximation, when 



Validity conditions for Eq. \5.21^ ; higher-order 
corrections 



(5.32) 



Note that this condition is equivalent to the requirement 
that, at least, one of the following two inequalities hold, 
either uj <^ £ or uj F . 

In the opposite case of large values of w, in the 
zero-order approximation, the term on the left-hand 
side of Eq. (j5.30|) can be neglected, yielding ry^"^ = 
a;~^arcsinh(2w/_F). Then, on expanding sinhwry in 
Eq. (I5.30p in a Taylor series in the vicinity of 77'^"^ and 
taking into account that [l^l 



arcsinha: = ln2a; — l/(4x^) -I- 0(a; "*), 
we obtain the expansion of 77 up to second order. 



(5.33) 



77 « - 



1 h ^ 1 In — 

uj uj^\ 2 F 



This series is rapidly converging when 



a;>£ln-, F. 
r 



(5.34) 
(5.35) 



All other parameters entering the present formulas are 
elementary functions of 77, and hence asymptotic expres- 
sions for them can be easily obtained using Eqs. (|5.3ip 
and (j5.34p . In particular, below we show asymptotic ex- 
pressions for 60 and Co defined in Eqs. (j5.28p and (|5.29l) . 
It appears that 60 increases with w. 



boK£ + F- 



FuJ^ 



2{£ + Fy 



2uj uj 
bQ ^ LO ~ £lii—-, a;>£ln— ,F; 

el' t 



(5.36a) 



(5.36b) 



whereas cq as a function of u increases for small uj and 
decreases for large uj, 



Fuj^ 



Co 



3(f -f F)4 

1 , 2uj 
Co « - In 



(5.37a) 



„- c^»£:in-,F. (5.37b) 
UJ et t 

Note that the asymptotic expansions for Ci are easily 
obtained from Eqs. (15.361) . since, according to Eq. (|5.28p . 
Ci = 5o - £. 

Below we will need a rough estimate of 5o for the whole 
range of the values of the input parameters £, F, and uj. 
The asymptotic expansions in Eqs. (j5.36l) imply such an 
estimate in the form 



Consider in detail the assumptions made in Sec. IV A II 
The above use of the stationary-phase method holds if the 
effective range of time integration in Eq. (|5.27p , which is 

on the order of 6q , is much less than tk+i —tk = tt/uj. 
In view of Eqs. (j5.36p . this requirement is satisfied under 
the condition 



w < 1. 



(5.39) 



In the usual units, this conditions means that uj <C 
\Eo\/h, i.e., the field frequency should be much less than 
the ionization potential divided by h. 

Furthermore, in Sec. IV A II we expanded the trigono- 
metric functions entering Eqs. (j5.4p and (15.51) in the small 
parameters, which is allowed if 



cj|4™'|<l {m>l). 



(5.40) 



Using Eqs. (15.171) and the limits for the values of the 
quantities in Eq. (|5.23p with the account of Eqs. ()5.3ip . 
(|5.34[) . (|5.36p . and (|5.37p . it can be shown that the in- 
equalities (|5.40p hold under the conditions (|4.13l) . (|4.14p . 
and ([Qg)) . 

To check the convergence of the series (|5.16p and (|5.19p 
for k = 21 {I = 0,1, . . .), we considered the asymptotic 
formulas for Eqs. ()5.17p and (j5.20p in the limits of low 
and high uj. We also obtained higher-order terms, up 
to fifth order, in expansions (|5.16p and (15.191) with the 
help of the software package Mathematica [3J] . The lat- 
ter expressions are rather unwieldy and are not shown 
here. The above analysis shows that under conditions 
(|n^, (liJi)) . and ((Og)) . the terms in the series (|ETC)) 
and (j5.19l) involving only z and q decrease rapidly with 
their order and, starting from third order, such terms can 
be neglected in (j5.27p . 

In contrast, terms involving pk require a special con- 
sideration. Indeed, Eqs. (|5.23p and (|5.37ap imply that 



\p2i\<{£ + Ff/{ujVF) (w«£ + F), 



(5.41) 



i.e., the width of the distribution for p2i increases with- 
out bound with decreasing uj, thus violating the above 
assumption of the smallness of p2/. However, for fixed uj, 
the terms proportional to a power of p2i decrease with 
the power. Correspondingly, our estimates of the terms 
in the expansion (|5.19p show that when uj is not too small. 



» {£ + Ff/F, 



(5.42) 



some third-order terms yield appreciable contributions 
and should be kept in the expansion, whereas the other 
third- and higher-order terms can be neglected. In par- 
ticular, we keep in the expansion (|5.19p (for k = 21) 
the third-order terms involving p2i, except for the term 
~ z'^P2i, which proves to be much less than the term 
~ 0^, so that we obtain 



bo--£ + F + ui. 



(5.38) 



fi,'Ht2i) ^ ~C2pIi/3 ~ q'p2l/{2bo), 



(5.43a) 
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where 



Co = 



g^(l + LoW) + 2gCi - 3gc^^?? + + 2^;^^ 
2^3 



(5.43b) 

Note that the frequency interval in Eq. (|5.42p coincides 
with the region, where the detachment rate depends sig- 
nificantly on uj [see the remark after Eq. (|6.1ip ]. 



B. Detachment rate and electron spectrum 

1. Momentum (Volkov- function) representation 
Inserting Eq. (|5.27p into (|5.3p . we obtain 

N 



fipit) = iy2^/^e'^''(*^'), (5.44) 

1=0 

where /p(i2;) = 7(^2; )U=o and 

/= [ df(t)oir)ze-''^-P+'-'">'^/'^. (5.45) 



An analysis shows that the interval — oo < z < <^o pro- 
vides a negligibly small contribution into the integral 
over z in Eq. (|5.45p . where is a number satisfying 

— 1/2 

1 <C Co ^ • (This result will be shown to be in- 

dependent of the specific value of Co-) In other words, 
/ is determined by large positive values of z, which 
reflects the fact that the electron tunnels in the posi- 
tive direction of the z axis when the force is positive, 
T{t2i) ^ £ + F > 0. 

Moreover, we assume, for simplicity, that the initial 
state is an s-state. This allows us to use the asymptotics 
()4.7p . which now can be further simplified, 

Mr) ~ {A/z) exphz - p^/{2z)] {z > 1). (5.46) 



Inserting Eq. (I5.46P into Eq. (|5.45l) yields 

= 2nA dzze-''«^'/2-g^-/2 

"'Co 
-Jl 



bo 



Here 



/i(e) = l-V^ee^'erfc(0, 



(5.47) 



(5.48) 



where erfc(^) is the error function [27]. To obtain the 
last equality in Eq. (|5.47l) . we replaced Co by 0, which 
practically does not change the result. 

Equation (|5.47p can be simplified as follows. As shown 
below, in Sec. IV B 41 the range of the values of for 
detached electrons is given hy q"^ ^ £ + F + to. There- 
fore, in Eq. ((CTT)) . the quantity q'^/VSki < {£ + F + 



u})b~^^'^ ~ {£ + F + u;)-i/2 < 1, where we took into 
account Eqs. (1535)) . and (01131). Thus, we ob- 

tain that q'^/^/Sbo ^ 1. Taking into account that 
/i(^) = 1 - + 2C^ -t- O(C^), we c an se t to a first 
approximation /i(g^/\/86o) « 1 in Eq. (|5.47p . yielding 



2TrA 

The detachment rate (12.41) becomes now 



W= lim 



dp 



(5.49) 



(5.50) 



where (3p{t) is given by Eq. (|5.44l) . The Volkov- function 
representation of the final state does not provide the en- 
ergy distribution of the detached electrons, which is di- 
rectly measurable in experiments. Moreover, this rep- 
resentation is not very convenient for calculations, since 
the summation in Eq. (|5.44p is difficult to perform ana- 
lytically, due to the imaginary terms proportional to t|; 
and t|j in the exponent in (I5.44[) [see Eqs. (|5.20p ]. There- 
fore, it is advantageous to transform the final state of the 
electron from the Volkov-function representation to the 
representation in the basis (I3.32p . as follows. 



2. Energy representation 



On inserting the function in Eq. (j3.34p written in the 
coordinate representation into Eq. (j5.1l) , the system wave 
function becomes 



^{r,t) = mMr)+ / dE 



dq 
(2^ 



■/3Eq{t)lpEqir,t). 

(5.51) 



Here 



dpz 
2-K 



hg{p,,E-Up)l3p{t). (5.52) 



On inserting here Eqs. (I3.26|) and (|5.44l) with the account 
of Eq. (|5.49l) and changing the variable of integration 
Pz p2i by Eq. (I5.15p . the terms proportional to t|; and 



t|; cancel out and we obtain 



2tt 



N 



X / dp2ie 



(=0 

.«C'3p|,/(6£)-coP2i/2-ia,p2i gg-j 



Here 



AE = E +1/2 = E - Eo 



(5.54) 



is the change of the average electron energy as a result 
of detachment and 



E 



diq' 



C3 = 1-2£C2, = ^ + (5-55) 
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where 



zi = —z- H — 



(5.56) 



The detachment rate of electrons with a given E, ir- 
respective of q, is obtained by integration over q in Eq. 
(f53T|) with the account of Eq. (f5J3l) . 



Equation (|5.53p can be further simphfied. We evaluate 
the integral in Eq. (j5.53|) with the help of the change 
of variable p2i p' = P2i + iEcq/C^, which eliminates 
the quadratic term in the integrand; then we sum the 
geometric progression. As a result, we obtain 



where 



2£_ 



1/3 



(27r)3/2 fdiq^/2-E' 



xe 



-C'5/2-E' /{2Ai)-C4,q^ /2+iTTAE/{2ui) 
2 _ ^2iiT{N+l)AE/uj 
I _ g2iTTAE/uj ' 



(5.57) 



where 



Here 



where 



dW 



Wn5{E-En), (5.668 



n— — 00 



Wn = 



Cs-K/Ai 



(5.66b) 



B„ = I Ai'f ) e-^du, (5.67) 



Ai 



C3 
2cn 



C5 — 2Co 



3cr 



A2 = 



1/3 



(5.58) 



(5.59) 



E'^E^E,, E,^Up-£zo, zo^zi-—^. (5.60) 
The detachment rate (|2.4p has now the form 



/•OO 



where the differential rate 



1 



dEdq t^oo t 



lim -IPeM . 



(5.61) 



(5.62) 



On inserting here Eq. (|5.57p . one can show that the rate 
of the detachment of an electron with given values of E 
and the transverse momentum q is 



dEdq 



= Wn{q)5{E~E^), 



(5.63a) 



n— — 00 



where 



■Wn{q) 



-Cz-E'^l Ax-Ctq^ 



xAi 



0<-'3 

A2 



(5.63b) 



= nw- 1/2 = n^iw + E^o, (5.64) 

i;; = -e,^ e„ -Up + Sz^. (5.65) 



The last expression in Eq. (j5.64l) is written in the usual 
units. 



d2 = di/{2Ci). 



(5.68) 



The total detachment probability is the integral of 
Eq. (|5.66ap over the energy, 



(5.69) 



n — — OO 



Equations ([537]) . ([5?63l) . (|5?66l) . and (|5?69| present the 
main results of the present paper. 



3. Asymptotic expressions 



In Sec. IV B 21 we introduced a number of new parame- 
ters. To get an idea about the values of these parameters, 
simple asymptotic expressions for them can be obtained 
with the help of Eqs. (jOTI) . ((534)) . (jOe)) . and ((07l) . As 
an example, we show for some of the above parameters 
the following simplified expressions valid in the limits of 
low and high frequencies. 



C3 



, F/2, uj<^£ + F 
^^"^^ oj/2, uj:$> £ln{uj/F), F 



F/{£ + F), uj<^£ + F 



1, 



uj > £ln(w/F), F 



(5.70) 



(5.71) 



F 1 



Flo^ 



-, a;<£ + F, (5.72a) 



a;2 2{£ + F) 8{£ + F)' 



1 ^ 

z„^--—, c^»£ln^,F. (5.72b) 



^^^^ uj-\ uj:$>£ln{uj/F),F ^^"^^^ 
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4- Discussion 



Let us discuss the above results. Equation (|5.63ap or 
(|5.66al) describes the spectrum of the detached electrons. 
These equations show that the energy distribution of the 
detached electrons consists of narrow equidistant peaks 
separated by hu. The energy En of the electrons in a 
given peak corresponds to absorption of n photons, if 
n > 0, or emission of —n photons, if n < [cf. Eq. ()5.64|) ] . 
The nth term in Eq. (|5.63ap or (|5.66ap provides the rate 
of increase per one atom of, respectively, the transverse- 
moment distribution or the total number of detached 
electrons with the energy £"„. 

The distribution of the electrons over the trans- 
verse momentum q is axially symmetric. According to 
Eq. (I5.63bp , the q distribution for a given energy is given 
by the Gaussian factor e~^'^'^ , which is independent of 
n, multiplied by the last factor in Eq. (j5.63bP . which is n- 
dependent. The factor e"'^*'^ shows that with the prob- 
ability close to one we have < C^^ ~ £ + F + lo, the 
latter relation being implied by Eq. (j5.73p . The asymp- 
totic expressions for the Airy function (I3.22p imply that 
the last factor in Eq. (j5.63bp can limit further the above 
interval for q^. Moreover, for positive E'^, this factor 
can provide oscillations of the number of the detached 
electrons as a function of q [cf. Eq. §2M]- 




FIG. 1. The spectrum, W„, in units of 2 x 10"^|£'ol/?i is 
shown by dots. Lines joining the dots are a guide for eyes. 
Here £ = 0.013, F — 0.05, and cj — 0.01 in the dimensionless 
units. 

The total energy spectrum integrated over 
Eqs. (|5.66p . also generally demonstrates oscillations re- 
sulting due to the static field, see Fig. [T] 



the quantity e , where Co is given in Eq. (|5.2ip with 
k = [and hence s = 1, as implied by Eq. (|5.26l) ]. In view 



of Eqs. ([5^ and (|5.20ap . the quantity 
the zero-order approximation to W, 



-2Co 



W (X e 



-2Co 



provides 



(6.1) 



This approximation is known as an approximation with 
"exponential accuracy" , since it involves a neglect of 
higher-order terms in the exponent in Eq. (j6.ip . It is 
interesting that Ec^. (|6.ip coincides with the formula for 
the rate of transmission through a one-dimensional tri- 
angular barrier in the presence of an ac field, obtained in 
a semiclassical approximation [61 . We consider the limits 
of Eq. (|6.ip for low and high frequencies, which will be 
used below. 

With the help of Eqs. KM . (lOTI) . Km . and (lOel) . 

we obtain the following simplified expressions. For low 
frequencies 



-2Co 



exp 



Fuj' 



3(£:^ 



F) 15{£ + FY 



3.2a) 



which agrees in the limit w — > with the exponential 
factor in the adiabatic expression ()4.12p and provides the 
first cj-dependent correction to the exponent. 

For sufficiently high frequencies, expanding Cq to the 
lowest orders in £ and F yields |3q . 



'2Co 



where 



-1/(2-) ^'££^ 

2uj J 



2M 



»i^,f2/3in^. (6.2b) 

r 



M 



1 

2a; 



4^3 



(6.3) 



In the last equality, we took into account Eqs. (I5.60p and 
(j5.72bp and the fact that Eq — —1/2 in the dimensionless 
units. Equation (|6.2bp hints that for high frequencies the 
detachment is caused mainly by niultiphoton absorption 
in the ac field [HI. (See Sec. I VI D] for more details.) 
We turn now to the exact results. 



B. The limit of a weak static field 



For a sufficiently weak static field, 
£<^F + uj, 



(6.4) 



VI. SPECIAL CASES 

A. The total rate with exponential accuracy 

As a preparation for the analysis of the above general 
formulas, we provide here the asymptotic expressions for 



the results obtained in Sec. IVB 21 can be simplified with 
the help of the following formulas. 



&o « 2d2 ~ l/Ci « Ca « 1, « 2Co, 

Aj^^ w (2/w) arcsinh(w/i^) - 2(F^ -f- u^Y^I^, 



13 



Using Eqs. (|6.2p and (|6.5|) . we obtain that in the hmit 
£ — >■ 0, expression (|5.63p reduces to the well known re- 
sult of detachment in an ac field [13, IsEl - ISTj , divided by 
two due to the present assumption (see Sec. IV A II) that 
we neglect the detachment occurring during the time in- 
tervals when the ac field is directed opposite to the static 
field. 



of Eq. (16.111) . Thus, this term provides an w-dependent 
correction to the adiabatic result (|4.12l) . Note that the 
condition for this term to be significant coincides with the 
validity condition (I5.42p of the present theory. Although 
for 



Fu^ <{£ + FY 



(6.12) 



C. The total rate for sufficiently low frequencies 

In this and following subsections, we provide simple 
analytical formulas for the total detachment rate in dif- 
ferent regimes. 

In the present subsection, we consider the case in which 
the spectrum ()5.66|) includes many peaks with compara- 
ble heights. As shown below, this occurs for sufficiently 
low frequencies, namely, when either of the following con- 
ditions holds. 



uj<^F or 



(6.6) 



In this case, the partial probabilities Wn slowly vary as 
a function of n, and the sum (|5.69p can be approximated 
by the integral 



UJ 



dEW{E), 



where 



W{E) = Wn\E,^^E- 



(6.7) 



(6.8) 



Using Eqs. ((6?8t . (j5.66bp . and ((5^ in ((6J)) and chang- 
ing the order of integration, we integrate first over the 
energy with the help of the formula. 



(a > 0), (6.9) 



derived in Appendix [BJ Performing the remaining inte- 
gration over u yields finally. 



W : 



2i]bo 



TT 
Co 



-2Co 



(6.10) 



Equation (|6.10p can be compared with Eq. (|6.ip . In con- 
trast to Eq. (|6.ip obtained with exponential accuracy, 
expression (j6.10l) provides the exact pre-exponential fac- 
tor. 

For low frequencies, Eq. (j6.10p can be simplified using 
Eqs. ([OT|) . (jOSaP . (jOTa)) . and (P?^ . 



W^^MiS + Ffe^p 



Flo' 



15{£ + Fy 



3(£ + F) 
uj<^£ + F. (6.11) 



This expression coincides with the adiabatic formula 
(|4.12p , if one disregards the second term in the exponent 



strictly speaking, the present theory is not valid, we can 
conclude from Eq. (j6.1ip by continuity considerations 
that in the region given in Eq. (I6.12p the total rate is 
described by the adiabatic result (|4.12p . at least, when 
uj<$:£ + F. 

Let us obtain the validity conditions for the above re- 
sults. The spectrum contains many peaks with compa- 
rable heights if the characteristic width of the spectrum 
Au, as a function of E'^ is much greater than the field 
frequency, 



3.13) 



Analysis of Eq. (j5.66bl) with the account of Eq. (|5.67p 
shows that, for sufficiently large and positive E!^, the 
spectrum is cut off due to the factor exp{—E!^/Ai), 
whereas for negative E'^ the dependence of Wn on i?^ 
is given roughly by 

Wn (X exp[|£;:j/Ai - (4/3)(K|/A2)--'/2], » A^. 

(6.14) 

As a result, we obtain that, for Ai ^ A2, the distribu- 
tion of Wn has a maximum near E'^ = and practically 
vanishes for E'^ ^ Ai and also for —E'^ ^ A2, so that 
the effective width of the distribution of Wn as a function 
of E'^ is A^ = Ai. Moreover, for Ai < A2, Eq. ((04)) 
implies that the distribution of Wn as a function of E'^ 
has a nearly Gaussian shape, 

Wn (X exp[-« - ^„,ax)VA^], Ai « A2. (6.15) 

This distribution is centered at the negative energy, 
-E'max ~ ^^2/(4^1), E^iid has the effective width A^ — 
Ag = ^y^JAl. The above results imply that for an 
arbitrary relation between Ai and A2 we have Ay; = 
max{Ai, Ag}. 

Consider two limits. For low frequencies, u) <ti £ + F, 
Eq. (|6l3l) is fulfilled. [Indeed, Eq. with the ac- 

count of Eqs. (|5.37al) and (|5.7ip implies that Ai = 
3{£ + FY /uP and hence we obtain that w ^ Ai < A^.] 
For high frequencies [a; > fhi(tj/F) , F] , Eq. ((535)) 
with the account of Eqs. ()5.37bp and (|5.7ip implies 
that Ai = cj/ln[2a;/(eF)], A2 = (£72)^^^, and hence 
Ag = {£Mn[2w/(ei^)]/(2tj)}i/2. Since now obviously 
a; > Ai, the condition (|6.13p is equivalent to a; ^ Ag, 
which can be recast as a; ^ ln(aj/F)]^/'^. Usually 
[ln(w/i^)]^/'^ ^ 1, and hence this factor can be omitted 
in the above inequality, yielding u) ^ E'^/^ . Thus, we ob- 
tain that Eq. (|6.13p holds under the condition that either 
Lo-^E + FoTLij<^ E'^l'^ ; this condition is equivalent to 
Eq. (EH). 
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D. High frequencies 



Consider the region of sufficiently high frequencies, 



> F, £2/3 In 



F 



(6.16) 



Now we can use Eq. (I6.2bp . Moreover, in view of 
Eqs. (15391) . (|5J0l) . and (jSJTt . now 



1/3 ^ 



(6.17) 



where is the characteristic frequency for the Franz- 
Keldysh effect [13-111; in the usual units [cf. Eqs. (|I3a|) 
and (|43c)) ] 



\2hmJ 



(6.18) 



In the present case, ^2 ^ A2 [cf. Eqs. (|6.16p and (16. 17^ ]. 
which allows us to simplify Eq. ()5.67p in two overlapping 
regions. 

First, for E'^ ^ A2, the square of the Airy function 
in the integrand in Eq. ()5.67p can be approximated for 
negative values of the argument by its asymptotic ex- 
pansion, the oscillations of which are smoothed out [cf. 
Eq. (|3.22bl) ]. and for positive values of the argument by 
zero [cf. Eq. p.22a|) ]. This yields 



Ai^ 



2fd2U-E'^\ E^'^e{E[^ - d.2u) 



(6.19) 



where 6'() is the step function. On inserting Eq. ()6.19p 
into (|5.67p and performing the integration, we obtain 
from Eqs. (|5.66bp . (glbl, (gSl), (ES]), and (IBTT]) that for 
sufficiently high-energy peaks, E'^ ^ ojg, i.e., n — M ^ 
ujg/uj, the partial detachments rates are 



2A^uj^/^F, (V2(n-M)) e-^/^^-) 



2lo 



where Fi(^) is Dawson's integral [27[, 



~ e ^ 



du •■ 



^«1 



1/(20, e»l 



(6.20) 



..21) 



The maximum of Fi{Ci is Fi(0.92) = 0.54. 

Second, for peaks with low and and moderately high 
energies, i?^ < or < i?^ ^ ^2, the exponential in the 
integrand in (|5.67p can be neglected, whereas the approx- 
imate expressions for the parameters used in the deriva- 
tion of Eq. (IS:^ [i.e., Eqs. ((OB . (|S3| and (PTT)) ] are 
still applicable now. As a result, Eq. (I5.66b[) becomes, 
for n - Af < 1 or n < M, 



-1/(2-) { ^\ 

2w J 



2n 



(6.22) 



where 



/oo 
Ai2(C)dC = ^{[Ai'(-0]' + eAi^(-0}- 

(6.23) 



Here Ai'(0 is the derivative of the Airy function 
The derivation of the second equality in Eq. (16.231) is 
given in Ref. [l3|. Using asymptotic expansions for the 
Airy function and its derivative \2f\ [cf. Eqs. p.22p ]. we 
obtain 



/(0«y! + 



cos(4g3/2/3) 
4e ■ 
exp(-4[e[3/V3) 



16[^[3/2 



-C » 1. 



i.24a) 
.24b) 



Moreover, Eq. (pT^ and the formula 10.4.5 in Ref. [27 
imply that 



/(O) 



32/3r2(i/3) 



0.210, 



3.25) 



where T{^) is the gamma function [27]. The first equality 
in Eq. ()6.23p implies that f{^) is a positive, monotoni- 
cally increasing function. 

In the region (j6.16p . detachment is a typical multi- 
photon process. The behavior is similar to that in the 
absence of the static field [17[ , but there are also impor- 



tant features peculiar for the present case of a non-zero 
static field. In particular, multiphoton detachment de- 
pends significantly on the threshold given by 

Eti, ^Mw^\Ea\+E,^ \Ea\ + Up~ £/u. (6.26) 

The threshold is increased by the value of the pondero- 
motive energy Up , Eq. (j3.17l) , which is a well known fea- 
ture of above-threshold ionization p7| . However, in the 
present case the threshold is also decreased by a quantity 
proportional to the static field. Moreover, in contrast to 
the case of a zero static field, now the threshold is not 
sharp, being smoothed by the quantity of the order of 
Wf. 

Equations (lOO)) and (p?^ together with Eqs. (|OT|) 
and (j6.24bp imply that Wn depends very strongly on n. 
The quantities Wn decrease rapidly with n increasing 
above the threshold, n — M > 1, and practically van- 
ish significantly below the threshold, M — n ^ Wf/w. 
However, Wn does not vanish for M — n ^ Wf/w, so that 
the effective threshold E'^.^^ is actually lower than Eth by 
an amount oi ^ ug. 

In a typical case, Wn is maximal for n — [Af]-|-1, which 
is the minimal number of photons required to excite the 
atom above the threshold. (Here [M] is the integer part of 
M.) Moreover, since the higher peaks decrease with n ap- 
proximately as a geometric progression with the common 
ratio equal to (eF)2/(2a;)2 < 1 [cf. Eq. (lOO)) ]. W[m]+i 
is much larger than the sum of the other quantities Wn, 
and hence the total detachment probability 



W !i^W, 



[A/J + l- 



(6.27) 



This relation can be incorrect only in relatively rare 
cases in which there is a peak sufficiently close to the 
threshold, so that the peak number no satisfies [no — 
Af [ <C 1 (where uq can be less or greater than Af). Then 
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the peaks no and Uq 
magnitude, so that 



1 can be of the same order of 



W w Wn 



(6.28) 



Figure [2ja) shows the partial rate Wq as a function 
of the frequency oj for two nonzero values of £ and in 
the absence of the static field. These plots show that the 
partial rate peaks when oj is near the threshold, i.e., when 
I no — A/ 1 <C 1. Therefore it is of interest to discuss the 
behavior near the threshold in more detail. In this case, 
detachment is described by Eq. ()6.22|) . which implies that 
detachment depends significantly on the static field. The 
static field results in several effects. First, the threshold 
energy is decreased in the presence of the static field [see 
Eq. (|6.26l) ]. Second, the threshold is smoothed out. As 
a result, detachment is possible below the threshold [see 
Eqs. (|6.22p and (j6.24bl) ] and is nonzero at the threshold 

[cf. Eq. m:m ]. 

Third, sufficiently above the threshold, the partial de- 
tachment probability, corresponding to the peak near 
the threshold, oscillates as a function of E'^^. Indeed, 
Eqs. and (|6.24ap imply that 



: COS 



3/2' 



L0£ < E'^^^ < LJ. 




0.080 0.085 0.090 
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FIG. 2. The partial rate, Wa, (a) and the total rate W (b) in 
units of 2 X 1Q~^^\Eo\/Tt, versus the frequency u) for F — 0.01. 
The solid lines: £ = 0.003, the dashed lines: £ = 0.001, the 
dotted lines: £ = Q. The quantities uj, F, and £ are given in 
the dimensionless units. 



The second term in Eq. (|6.29p demonstrates oscillations 
oiWno as a function £'^^ = a;(no — M). The quantity £'^^ 
can be changed, and hence the oscillations of Wng can be 
observed on varying any of the quantities cj, F, and £. 
The oscillations decrease with i?^^ ; the amplitude of the 
oscillations and the intervals between them increase with 
£. Note that Wn^ decreases with to sufficiently above the 
threshold [cf. Fig. ©(a)]; this feature is determined by 
the factor before the braces in Eq. (|6.29l) , while all other 
frequency-dependent factors in Eq. (|6.29p increase with 

LU. 

Since the partial rate Wno can provide a significant or 
even dominant contribution to the total rate (16.281) . the 
above features of the near-threshold behavior generally 
apply not only to Wno but also to W, see Fig. ©(b). 
The above near-threshold behavior is reminiscent of that 



obtained in the case of the Franz-Keldysh effect [401 - 43 
or one-photon detachment of ions in an electric field 11, 
[T3|. There are, however, also features peculiar to the 
present case. In particular, in contrast to the Franz- 
Keldysh effect, in the present case, oscillations of the 
detachment rate can occur, as mentioned above, not only 
as a function of the laser frequency, but also as a function 
of the laser-field intensity or the static-field strength. 



VII. APPLICATION TO ELECTRON EMISSION 

Consider electron emission from a metal or semicon- 
ductor surface or from photosynthetic bio-complexes due 
to a static field perpendicular to the surface, under con- 
ditions typical for scanning tunneling microscopy 043- 
Let Vb be the tip (positive) potential with respect to the 
surface and d be the tip-surface separation. Then the 
static-field strength £f ~ Vo/d. In view of Eq. ()4.3cp . 
the dimensionless quantity £ = \eo\£f/Fa or, in view of 
Eq. gjg). 



^ 0.0974 



d\Eo\^/' 



(7.1) 



where d, Eq, and Vq are measured in nm, eV and V, 
respectively. For example, for d = 1 nm, Eq = —5 eV, 
and Vb = 1 V Eq. ^ yields £ = 0.0087. 

Assume that a laser (or ac) field perpendicular to the 
surface is superimposed on the static field. As follows 
from the above theory, the laser field affects electron 
emission only if the field intensity is sufficiently high. 
Let us estimate the minimal strength of the laser field re- 
quired to affect significantly electron emission. The lower 
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limit on the ac-field amplitude, at least, for low frequen- 
cies, Lo <t: F + £,is implied by Eq. (HJ5)) . F > f^, or 

Fu>£a£^, (7.2) 

where Fu is the laser electric field amplitude in V/cm, 
the subscript u indicates that the quantity is measured 
in usual units, and 

£a = i^a/|eo| = 1.03 X 10«|So|3/2 v/cm. (7.3) 

Here Eq is in units of eV, and in the last equality we 
used Eq. (|4.3cp . For example, for Eq = —5 eV, we obtain 
£a = 1.1 X 10^ V/cm. If also 8 = 0.01, Eq. jLl]) yields 
Fu ^ 10^ V/cm, which corresponds to the laser inten- 
sity II = fscF'^j2 > 1.7 X 10'^ W/cm^, where eo is the 
vacuum permittivity. Note that the minimum laser in- 
tensity in the latter condition increases with \Eq \ as |i?oP 
[cf. Eqs. dZSl) and iTlh ]. 

Consider the field frequency. The present theory im- 
plies that electron-emission rate increases with the laser- 
field frequency when the latter is sufficiently high. Let us 
estimate the minimal frequency required to affect signif- 
icantly electron emission. Equation (I4.3ap implies that 
the field frequency measured in hertz, w„, is related to 
the dimensionless frequency w by 

a;„ = 4.84 X lO^l^olw Hz. (7.4) 

As follows from Eq. (|6.1ip . the field frequency signifi- 
cantly affects the emission rate for 

. > ^il (7^5) 

Vf 

[cf. also Eq. ^5A2\ ]. If £ = F = 0.01 and Eq = -5 
eV, then Eqs. (j7.5p and (|7.4p imply that the frequency 
affects the transition rate for lUu > 10^^ Hz. In view of 
the validity condition (|5.39p . we obtain in this example 
that in the frame of the present theory, the frequency 
affects detachment in the interval 10^^ < < 2 x lO^'' 
Hz or 1.5 ^ A < 30 /im, where A is the wavelength of the 
laser field. 



VIII. CONCLUSION 

In this paper we have studied multiphoton electron 
detachment from atoms or negative ions by a sum of a 
static and a laser (or, more generally, ac) fields. Though, 
strictly speaking, the present theory holds for electron 
detachment from negative ions, the main qualitative fea- 
tures of detachment revealed here should be applicable 
also to ionization of neutral atoms. We have developed an 
analytical theory of this phenomenon and derived simple 
formulas valid in a broad range of the parameters of the 
problem. In particular, we have obtained that the energy 
spectrum of the electrons consists of narrow equidistant 
peaks separated by huj. 



We have identified, at least, two physically different 
regimes. For low and moderately high frequencies, there 
are many spectral peaks of comparable size, and, as a 
result, the total rate smoothly depends on the frequency. 
For very high frequencies, the detachment is reminiscent 
of that in the absence of the static field, including the 
appearance of the energy threshold. We have shown 
that the static field can significantly affect detachment 
for both low and high frequencies. The effects of the 
static field include, in particular, oscillations of the en- 
ergy spectrum of the products for low frequencies and the 
decrease and the smoothing of the threshold and oscilla- 
tions of the total rate for high frequencies. The results 
of the present work significantly clarify the physical pic- 
ture of multiphoton electron detachment in the presence 
of a static electric field and provide a basis for further 
investigations in this field. 
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Appendix A: Average energy in different states of 
an electron in static and ac fields 

Here we calculate the quantum mechanical expectation 
values of the energy in different states describing a free 
electron simultaneously affected by static and ac fields 
[Eq. (|3.2p ]. More specifically, we shall consider the Volkov 
functions and the functions \4>Eq{t))- We shall also obtain 
time averages of the expectation values of the energy over 
the electron oscillations. 

The energy operator of an electron in a static field is 
given by [in the units p.3p ] 

H=—-8-r = K + U, (Al) 

where K and U are the kinetic and potential energies, 
respectively. An ac field can change the electron energy 
with time. It is of interest to obtain the instantaneous 
and time-averaged expectation values of the electron en- 
ergy in the presence of the ac field, F sin wt, for the basis 
states that are of relevance here. The expectation value 
of the electron energy is given by the average of the op- 
erator H [Eq. (lAip ] in the electron state. It is convenient 
to work in the momentum representation, in which the 
kinetic and potential energies of an electron are, respec- 
tively, 

P2 ^ d 

K^—, U^-iE-—. (A2) 
2 ' dP 
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A difBculty in obtaining expectation values for basis 
states of the continuum is that such states are not nor- 
malizable to one [cf. Eqs. ()3.7|) and p.31|) ]. To overcome 
this difficulty, we consider expectation values for wave 
packets and obtain the result for a basis function in the 
limit when the width of the wave packet tends to zero. 



The expectation value of the energy in the state (|A9 



IS 



-Ecxp — {H) 



dP - 
-^^Sl{P,t)H^,{P,t) = {K) + (U), 

(AlO) 



where, in view of Eqs. (jA2[) and (jA9l) . we obtain 



1. Volkov functions 

We begin with the Volkov functions. The Volkov func- 
tion in the momentum representation. 



(A3) 



is given in the general case, in view of Eq. (j3.4p . by 

^p„iP,t) = (2^)35(P- Po(t))e^/o (A4) 
Here 



P^{t)=P^ 



eo 



Ait), 



(A5) 



where A{t) is given in Eq. p.9p . It is easy to see that the 
functions (jA4[) satisfy the orthonormality relation in the 
momentum representation [cf. Eq. (13. 7|) ]. 



{P') 3Ag P^jt) 
2 2 2 ' 



(Alia) 



([/) = 



-£■ dT 
Jo 

-£■ f drPiir) 
Jo 



P) + -Ait) - -Air) 
c c 



(Allb) 



Here in the last equality we took into account the equality 
(P) = Piit) and Eq. (jM]) . 

The expectation value of the energy for the Volkov 
function is obtained in the limit Ap — > 0. Then, on drop- 



ping the subscript 1 and taking into account Eqs. (|3.5|) 
and p.9|) . we obtain from Eqs. (|Alip that for the Volkov 
function \ippit)) the mean kinetic and potential energies 



P^it) 1 / ^ - F ^ 
{K) = = 7:\P + £i cos wt 



dP 



3 iP, OV-p- iP. t) = i2iT)Hipo - p'o)- (A6) 



We consider a Gaussian wave packet in the momentum 
space. 



f27r)~^/'' 
*i(P,0 = ^ / dpocxp 



iPo -Pif 



4A2 



i^PoiP.i)- 
(A7) 

Here pi is the average value of the vector po , whereas Ap 
is the standard deviation of the components of pb- In the 
limit Ap — >■ 0, the function in Eq. (jA7l) becomes propor- 
tional to the Volkov function ijjp^ (P, t) . It is easy to see 
that, in view of Eq. (|A6[) . the state (jA7p is normalized to 
one, 



-p-£t- 



£ ■ F 



■ sin ujt. 



(A12) 



Equation (|A12p shows that the magnitudes of the mean 
kinetic and potential energies increase as t^. However, 
in calculating the expectation value of the total energy 
Ecxp — (K) + ([/), the terms proportional to and t 
cancel, and we obtain 



2 



p- F 



cos Lot 



F' 



^£-F 



sin ujt 



2^2 
t cos iot 

UJ 



cos^ Ujt 



(A13) 



Finally, averaging E^xp over the electron oscillations, we 
obtain the simple result. 



dP 



\^i{P,t)? 



(A8) 



On inserting Eq. (KM into Eq. (K7\ . we obtain *i(P, t) 
in an explicit form. 



dr 



4A2 



P+-Ait)~-AiT)y}, (A9) 
c c J 



where Pi it) is given by Eq. (jA5|) . 



E, 



cxp 



EL 

2 



(AM) 



Consider two important special cases. First, for the 
familiar Volkov functions describing an electron in a si- 
nusoidal field, the above formulas hold with £ = and 
U = 0. As a result, Eq. (|A12p implies that 



E, 



1 



cxp = = ^ I coswi 



(A15) 



The time-averaged expectation value of the energy is 
again given by Eq. (IA14p . 
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Second, for the Volkov functions describing an electron 
in a static uniform field, we have F = 0. Then Eqs. (|A12[) 
and (|A13I) yield 



(K) 



(U) 



~p-£t 



2j-2 



2 ' 

(A16) 



2. Functions lipEqit)) 

We shall consider the functions \ipEq(t)) for the case 
of interest here — when the ac field is linearly polarized 
along the static field, i.e., when Eqs. p.lOp and p. lip 
hold. The representation of these functions in momentum 
space is obtained as in Eq. (IA3p . So, in view of Eq. p. 301) , 
we obtain 



^hl{po,,E~Up)i^p„{P,t), (A17) 



where po = {pox,Poy,Poz) and 9*0 = {Pox,Poy)- Inserting 
Eqs. (IX^ and (jM)) [where Po{t) is defined in Eq. (|X5|) 
with A{t) given in Eqs. p.l2p and (j3.13l) ] into Eq. (|A17p 
and performing the integration yields 



exp • 









6 



dTP^{t,T) 







d{Q - qo), (A18) 



where E = E - Up ~ P,{t) = P, + (eo/c)^(t), 

P,{t,T) = P,{t) - (eo/c)A(r), and Q = {Px,Py). The 
function (jA18P can be obtained also by a direct solution 
of the Schrodinger equation [§]. 

Consider a wave packet in the momentum space with 
a Gaussian distribution of the quantities E and qo, 



and performing the integration yields 

25/4^3/4^1/2 C A|F2(^) 



■ exp ■ 



4A2 



E 



£2 



Ei-U-o 



2 



iQH 



pm 

6 

(A20) 



Let us obtain the expectation values of the kinetic and 
potential energy in the state (IA20p . Taking into account 
the expressions for the kinetic and potential energy given 
in Eq. (jA2p . where now U has the form —i£-^, we ob- 



tain 



{K) 



8A| 



2c2 

2 



{U)=E,-Up-^l-^- 



,2 + i 
£2 



8A| 



dTA(r). 



(A21) 
0, since 



Both {K) and {U) diverge in the limit A^ 
they contain terms ~ A^^. 

However, the expectation value of the total energy 
-Eoxp = {K) + ([/) is independent of both A_b and Ap. 
Dropping the subscript 1, we obtain that the expectation 
value of the energy in the state \'4'Eq{t)) is 



EcKp — E Up 



I dTA{T) (A22) 



sincLit t cos Lot 



or, taking into account Eq. p.l3p . 

F2 ^ ^ 

Ec^r, = £■ - f/„ + — r cos^ ut + E-F 

2w^ \ uj^ uj 

(A23) 

Finally, the average of this quantity over the electron 
oscillations yields 



Ecxp — E. 



(A24) 



(2^) 



-7/4 



(go - gl)' 
4A2 



dE / dqo exp 



4A| 



(A19) 



Here £^1 and (fi are the average values of E and qo, re- 
spectively, and Ae (Ap) is the standard deviation of 
E (of each component of go)- In the limit Ae,Ap ~^ 
0, the function in Eq. (jA19P becomes proportional to 
i^Eiqi{P,t). Using Eq. p.3ip . it is easy to show that 
the wave function in Eq. (|A19p satisfies the normaliza- 
tion condition given by Eq. ()A8P where ^i{P,t) is re- 
placed by *2(^,0- Inserting Eq. (|JT8)) into Eq. (UT9)) 



Appendix B: Calculation of the integral in Eq. (|6.9p 

To calculate the integral in Eq. (|6.9j . we use the dif- 
ferential equation for w{^) = Ai^(^) [23|, 



w 



Let us denote 

J = 



A^w -2w^0. (Bl) 



d^Ai^{^)e''^ (a>0). (B2) 



Multiplying the both sides of Eq. (IBip by and inte- 
grating each term over ^, using the rule of integration by 
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parts for the first two terms, we obtain the differential 
equation for J, 



'da ~ VT ^ 2a) 
Equation (jB3p can be easily solved, yielding 



J 



(B3) 



(B4) 



where C is an unknown constant. 

To obtain C, we note that Eq. (|B4p implies that J 
diverges in the limit a — > 0, so that 



(B5) 



The asymptotic equations p. 221) imply that the divergent 
term in Eq. (IB2p arises due to the integration over large 
negative values of i^, i.e., for a — > [cf. Eq. (j3.22bp ] 



2vrViei 



l-sin(||er 



(B6) 



where is any number such that ^ 1- The replace- 
ment of ^0 by in Eq. (jB6| changes J by a term which 
is finite for a — >■ 0. Moreover, the integration of the oscil- 
lating term in Eq. (IB6p can be shown to yield the result 
which is finite for a — > 0. Terms, which are finite for 
a — ?> 0, are of a higher order than the terms ^ a^^/^ and 
hence can be neglected to a first approximation. There- 
fore, in the limit a — >■ 0, Eq. (|B6p implies that 



J-oo 2Try/\^\ TT Jo 



(B7) 



where the first integral was transformed by the change of 
variable ^ = — C^. A comparison of Eqs. (jB5|) and (jB7[) 
yields the equality G = (2-y/7r)~^, inserting of which into 
Eq. ((B4)) results in the formula 
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